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ABSTRACT 

The  fluctuation-dissipation  formula  for  the  viscosity  in  terms 
of  the  stress  tensor  fluctuations  reproduces  the  experimental 
logarithmic  temperature  dependence  of  the  hydrodynamic  viscosity. 
Using  the  theoretical  wave  number  and  frequency  dependent  viscosity 
in  the  problem  of  critical  diffusion,  we  find  that  the  effects  of 
non-locality  and  retardation  practically  cancel,  resulting  in 
satisfactory  agreement  with  a  "•ecent  light  scattering  determination 
of  the  effective  viscosity. 
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Kawasaki1  has  shown  that  the  critical  variation  of  the 
diffusion  coefficient  in  a  binary  liquid  near  its  critical  point 
can  be  represented  by  the  Einstein  relation  D  =  T  u^j  where  T  is 
the  temperature  (we  use  natural  units  such  that  Boltzmann's  constant 
is  uni*y) ,  £  is  the  con-elation  length  for  the  concentration 
fluctuations,  and  -  (6i-ri£)  1  is  Stokes'  formula  for  the  mobility 
of  a  sphere  of  radius  £  moving  through  a  liquid  of  viscosity  n . 

2 

This  result  has  also  been  established  by  one  of  the  present  authors 
by  a  different  method,  and  gives  the  rate  of  relaxation  of  a  con¬ 
centration  fluctuation  of  wave  number  q  as 


VDq 


2  =  Tq« 

6rrn£ 


(1) 


provided  £<<.q-1;  i.e. ,  the  wavelength  should  be  much  greater  than 
the  correlation  length.  But  as  the  critical  point  is  approached, 
T+’i  and  £+°°.  The  above  inequality  is  then  no  longer  satisfied  and 
D  becomes  a  function  of  q,  corresponding  to  "non-local"  diffusion. 
This  change  is  carried  out  in  Eq.  (1)  by  substitution  of  an 
effective  value  for  £_1  according  to 


^  *eff  "  aef&* 

where  a^^  is  some  numerical  constant  of  order  of  magnitude  unity. 

Eq.  (2)  expresses  in  a  quantitative  form  the  qualitative  rule  of 
3  4 

dynamical  scaling  ’  that  all  temperature  dependence  is  to  be 
expressed  in  terns  of  £  and  that  as  T-VT  all  factors  of  £  become 
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replaced  by  the  wavelength.  Actual  computation  yields 

a  cc  =  3tt/8  so  that  in  the  limit  T-+T  we  obtain 
eff  c 


c  Tq 

Fq  =  "l6neff(qT 


(3) 


5  6 

In  Eq.  (3)  we  have  allowed  for  the  fact  that  recent  experiments  ’ 
have  clearly  established  a  critical  temperature  dependence  in  the 
hydrodynamic  viscosity  n.  It,  therefore,  as  well  as  the  explicit 
factor  ?  \  has  to  assume  a  q-dependent  effective  value,  which  we 
denote  by  n^^Cq).  The  purpose  of  the  present  note  is  to  compare 
our  calculation  of  this  function  with  some  measurements  of  r 

.  q 
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recently  carried  out  in  this  laboratory.  As  the  theory  is  rather 

involved,  we  prefer  to  present  our  results  in  a  semiphenomeno logical 

fashion,  touching  only  the  main  points  of  the  calculation  at  the  end. 

8 

Proceeding  in  this  spirit  we  first  take  note  of  the  separation  of  the 
hydrodynamical  viscosity  into  an  "ideal,"  non-critical  background 


n.^  and  a  residual  critical  portion  An.  An  accurate  empirical  fit 
to  the  ratio  of  the  critical  portion  to  the  total  as  a  function  of 
5  for  the  binary  liquid  3-methylpentane-nitroethane  is 
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An(g  L) 

nCT1) 


=  A  In  (qD0> 


(4) 


where  A  =  0.051.  The  correlation  length  is  known  from  light  scattering 
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intensity  measurements  to  have  the  temperature  dependence 


(£;/£  )  =  (T-T  /'T  )  v  with  C  =  2.56A  and  v  =  0.616.  New  sub- 
o  c  c  o 

vise 

stituting  from  Eq.  (2)  and  denoting  the  numerical  factor  by  , 

we  obtain  at  the  critical  point 


Wq^ 

neff(q) 


id 


=  A(ln 


qD 


i  vise  v 
-  In  a  «  ) 
eff 


(S) 


The  fact  that  the  functional  dependence  is  logarithmic  permits  us  to 
consider  the  case  that  the  logarithm  is  very  large.  The  constant  term 
can  then  be  neglected.  To  this  "logarithmic  accuracy,"  we  expect 
the  fractional  critical  viscosity  which  is  effective  in  the  diffusion 
process  to  be  precisely  the  same  function  of  q  that  the  hydrodynamic 
viscosity  is  of  This  function  is  plotted  as  the  solid  line  in 

Fig.  1.  The  values  of  neff(q)  determined  experimentally  from  line 

widths  are  shown  as  circles.  The  parallel  dashed  line  illustrates 

•  •  vise  ’-I  ♦ 

better- than- logarithmic  accuracy  for  ag^ '  =  2  ,  which  would  follow 

from  an  incomplete  theoretical  treatment,  as  explained  below.  This 

line  serves  to  indicate  the  sensitivity  of  the  logarithmic 
.  •  v^  sc  .  . 

approximation  to  ag^  .  It  is  evident  in  Fig.  1  that  for  the  theory 

♦  vise 

to  agree  with  the  data  the  computed  value  of  ag^  should  no.-  be 

significantly  smaller  than  unity.  The  remainder  of  this  note  is 

devoted  to  establishing  that  this  is,  indeed,  the  outcome.  In  fact, 

we  find  from  a  complete  theoretical  treatment  (i.e. ,  one  taking  into 

account  both  the  wave  number  and  the  frequency  dependence  of  the 

viscosity)  =  0.92.  This  is  so  close  to  unity  that  there  is  no 


L 


need  to  draw  an  additional  line.  We  can  regard  the  hydrodynamic 

function  (a  =  1.00)  as  a  sufficiently  accurate  representation  of 

the  theoretical  effective  viscosity  and  we  see  from  Fig.  1  that  it 

vise 

is  consistent  with  the  experimental  finding  ag^  (EXP)  =  2±1. 

In  order  to  achieve  better-than-logarithmic  accuracy-,  we  need 

•  •  •  •  vise 

a  theory  of  critical  viscosity  for  calculating  ae£f  *  Such  a 

.9  *  .10 

theory  has  been  proposed  by  Kawasaki  and  by  Deutch  and  Zwanzig 

We  have  followed  the  latter  approach,  which  is  based  on  the 

11  12 

fluctuation-dissipation  theorem  ’  connecting  the  viscosity  with 

fluctuations  in  the  off-diagonal  components  of  the  stress  tensor. 
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In  this  way  we  have  obtained  an  expression  of  precisely  the  form 

.  14 

of  Eq.  (4)  with  a  theoretical  value  for  the  coefficient  A  equal 

2 

to  8/15it  =  0.054,  in  close  agreement  with  the  experimental  value 

0.051.  (The  Debye  cutoff  q^  is  a  free  parameter  in  the  theory  and 
is  to  be  fixed  by  taking  it  equal  to  the  experimental  value. )  The 
fluctuation-dissipation  theorem  expresses  the  wave  number  and 
frequency-dependent  viscosity  as  a  four-dimensional  space-time  Fourier 
integral  over  the  correlation  function  of  fluctuations  in  the  off- 
diagonal  component  T^  of  the  stress  tensor,  (x  and  y  are  particular 
Cartesian  coordinates . )  It  is  more  convenient ,  however ,  to  treat 
the  time  dependence  by  the  Laplace  transform,  in  which  we  write  the 
"frequency"  as  times  the  dimensionless  variable  y  .  Thus  the 
generalized  viscosity  is 


i(q,y) 


2T 


J 


2 

4  q'*21  ~  1^21^ 

d  21  e  q  <T  (?)  T  (1)>,  (6) 

xy  xy  ’ 
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where  the  integration  is  over  all  of  space-time,  q  is  in  the  y- 
direction,  and  the  angular  brackets  denote  the  thermal  equilibrium 
average.  We  are  interested  in  the  critical  fluctuations  in  T  , 
which  are  given  by  the  canonical  expression3-3 


x  3s 


s  s  , 
x  y 


where  the  Landau-Ginzburg  free  energy  density  depends  quadratically 
on  the  partial  derivations  of  the  concentration  s^  through  the 
gradient  term  (Vs  )2/ 22.  Inserting  Eq.  (7)  into  Eq.  (6)  and  factoring 
the  ensemble  average  in  terms  of  the  concentration  correlation 
function  * 


G(21)  =  <  s(2)  s(l)  >  =  (2tt)_3  jd3q  e 


-3  fj3  J-  S‘*21  “  rq|t21l 


g(q)  (3) 


enables  us  to  carry  out  tie  integration  in  terms  of  the  Fourier 


transform 


g(q)  = 


Substituting  Eq.  (7)  into  Eq.  (6)  and  factoring  the  fluctuating 
concentration  variables  into  pairs  we  find  for  the  critical  viscosity, 
after  a  certain  aiTount  of  manipulation, 

64ir  q5 


X 
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2  2  ^ 

g(q1)g(q2)^1‘q2> 


d,ldq2<ll<>2  rc  tfc  trcy3 


2,  2^  2 


2  2, 


[2q  (q1+q2)-q  -(q1~q2>  3 


q1+q2^q 

Iq^qjUq 

ql,2~qc 


Tq° 

647i2rc 

q 


q2  q 


du^duj 

V2 


Ul+U2"l 

luj-UjUl 

“l.iV'l 


2  2  ^ 

(uru„)  2  7  o  o  2 

"X  3 .  3  ^(u^Uji-l-Cu^Uj)  ] 

u^+u2+y 


(10) 


The  cutoff  qc  is  proportional  to  the  Debye  cutoff  q^.  The  static 
case  y  =  0  can  be  integrated  analytically  and  yields 


An(q ,o)  _  _ 8_ 


'eff 


(q) 


15tt^ 


(In  4  13  4tL_) 

q  15 


^  (In  -—+  0.062).  (11) 

9/1  IBtt2  q 


Eq.  (10)  has  been  written  for  the  special  case  £  ^=0.  Finite  values 

-2  -1  -21/2 

of  £  require  correction  factors  [l+(q£)  ]  ard  [l+(q£)  J  in 


g(q)  and  r  ,  respectively 
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Comparison  of  the  result  for  the 


q  =  0,  E,  1>0  hydrodynamic  case  with  the  q>0,  E,  ^=0  non-local  case 
gives,  in  the  language  of  Eq.  (2) 


N.L. 

Wf 


2  exp  (- 


4tt 
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9/1 


)  =  0.496. 


(12) 


If  we  were  to  assume  that  the  frequency  dependence  of  the  viscosity 
could  be  neglected  in  the  critical  diffusion  process,  we  would  have 


8 


. ,  ...  vise  - 

the  approx  urat ion  .  a 

by  the  dashed  line  labelled 


N.L. 
eff  ‘ 

"N.L. 


II 


This  approximation  is  depicted 
in  Fig.  1.  It  is  interesting 


17 

to  note  that  Eq.  (12)  confirms  the  general  rule  of  thumb  that 
when  the  factorization  brings  in  the  correlation  length  twice  (via 


the  equal-time  Green's  function),  ae^;l/2. 

Biit  Eq.  (12)  cannot  directly  be  applied  to  the  diffusion 


problem  because  it  is  applicable  only  in  the  static  limit.  The 
approximation  a^^.c-  a^^'  is  too  crude,  as  the  diffusion  process 
depends  'jpon  the  relaxation  of  current  fluctuations,  which  are  the 


product  of  concentration  and  velocity  fluctuations.  Because  of 


the  finite  concentration  relaxation  rate  we  are  forced  to  study 


viscuous  damping  of  the  velocity  field  at  non- zero  frequencies. 

Ihis  results  in  a  kind  of  "retardation"  correction  to  the  critical. 

viscosity  which  leads  to  a  decrease  in  the  effective  viscosity. 

Our  final  task  is  the  computation  of  this  decrease,  which  entails 

3 

the  use  of  Eq.  (10)  at  all  frequencies  y  .  We  define  a  dynamical 
scaling  function  o(y)  by 


n(q,Y)~ri(q,o)  _  _ §_ 

neff (q)  15tt2 

For  0<y<<l,a(y)  =  -0.316y^ 
approximation 


o(y)*  (13) 

,  while  for  y>>l  we  have  the  high  frequency 


o(y )  =  In 


,1/ 3 


tt/3  ,02/3  s-2  13  .  2w/3 _ 

sunCn/’S  ^  Y  15  3sin(iT/3) 


0.169  -  Iny  -•  0.479Y 


-2 


(14) 


We  have  confirmed  by  numerical  computation  that  the  above  low 
and  high  frequency  approximations  are,  in  fact,  very  accurate 
for  y<0.5  and  y>l. 5 ,  respectively.  Interpolation  through  the 
region  y=l  yields  a  smooth  negative-definite  monotonic  curve. 
Therefore,  the  effect  of  retardation  is  necessarily  to  decrease 
the  logarithmic  divergence  by  a  weighted  mean^,  o^,^-,<0.  To 
determine  c^ve  we  exPar'ded  the  reciprocal  of  the  viscosity,  which 
appears  in  the  time-dependent  velocity  correlation  function,  in 
powers  of  the  critical  portion  An,  and  worked  to  first  order  in  An* 

The  resulting  two-dimensional  integral  over  the  wave  numbers  of  the 
concentration  and  velocity  fluctuations  can  be  reduced  analytically 
to  quadrature  in  terms  of  their  ratio  y.  Using  'the  values  of  o(y) 
obtained  as  described  above,  and  carrying  out  this  integration 
numerically,  we  have  found  0^^  -  -0.62,  corresponding  to 

vise  N.L.  °AVE  n  n~  ,,, 

a  *:*:  =  a  rr  e  -  0.92.  (15 

eff  eff 


As  mentioned  above,  this  is  in  satisfactory  agreement  with  the 
experimental  value  a^-^C  (EXP)  =  2±1.  The  fact  that  this  coefficient 
comes  out  so  close  to  unity  is  a  consequence  of  two  competing  phvsical 
effects,  namely  the  non- locality  and  the  retardation. 

A  further  outcome  of  the  theory,  which  we  mention  onlv  in  passing 
is  that  the  retardation  produces  some  non-Lorentzian  distortion^  in 
the  diffusion  line  shapes  which  may,  however,  be  too  weak  an  effect  to 
detect  experimentally.  Finally,  it  is  a  pleasure  to  acknowledge  a 
stimulating  conversation  with  Professor  L.  Kadanoff  and  many  helpful 
discussions  with  R.F.  Chang,  P.H.  Keyes,  J.V.  Sengers,  and  C.O.  Aliev. 
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CAPTIONS  FOR  THE  FIGURES 


Fig.  1  Fractional  critical  viscosity  vs.  wave  number.  The  solid 

line  labelled  "HYD!1  represents  the  experimental  logarithmic 
dependence  of  the  hydrodynamic  viscosity  upon  the  inverse 
correlation  length  (references  5  and  6).  The  dashed  line 
labelled  "N.L."  shows  the  static  non-local  correction  for 
the  theoretically  computed  value  of  a^  L  =0.496.  This 
would  be  the  result  for  a  crude  theory  which  did  not  include 
retardation.  In  the  complete  theory  the  correction  for 
retardation  cancels  the  effect  of  static  non-locality  and 
brings  the  theory  back  down  to  the  solid  line.  Thus  the 
solid  line  represents  both  the  hydrodynamic  viscosity  and 
the  theoretically  expected  effective  viscosity.  This  line 
is  in  satisfactory  agreement  with  the  circles,  which  show 
the  light  scattering  line  widths  measured  by  Chang,  Keyes, 
Sengers,  and  Alley  (reference?). 


